Grey solitons in doped nonlinear fibers 
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Grey solitons, with complex envelope, are obtained as exact solutions of the coupled system 
describing doped optical nonlinear fibers. The phase degree of freedom plays a crucial role in 
removing the restrictive equal frequency conditions on the tanh-sech paired pulses, identified earlier 
0|. The cross-phase modulation is found to control the speed of the grey soliton, which has an upper 
bound. The coupled soliton solution obtained here was found to be stable in a wide parameter range. 
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I. INTRODUCTION 

Self-induced transparency is a dramatic manifestation 
of nonlinear effects, leading to coherent localized pulse 
propagation in an atomic medium The controlled, 

localized excitation and de-excitation of an ensemble of 
two level systems is effectively captured by the hyperbolic 
secant soliton pulse profile, which propagates as a stable 
soliton. The Maxwell-Bloch equations aptly describe the 
pulse and averaged atomic population dynamics, which 
have been extensively investigated for multi-level systems 
SIMS- Localized and nonlinear cnoidal wave solutions 
have been identified, as exact solutions [f| 0] , which has 
found experimental verification [Til [l2| . 

The subject of pulse propagation in nonlinear fibers 
has also been well-studied, since its prediction by 
Hasegawa in 1974 [l3j]. The non-linear Schrodinger equa- 
tion (NLSE), where the cubic nonlinearity arises due to 
Kerr effect, describes the pulse dynamics. This integrable 
system in one dimension possesses soliton solutions 
which can be dark [l|| , bright [l7j - r20j or grey 



depending on the nature of nonlinearity and dispersion. 
Some of these solitons were experimentally observed, long 
after their prediction (22[. Gross-Pitaevskii equation in 
one-dimension, which captures the mean-field dynamics 
of Bose Einstein Condensates (BEC), being identical to 
NLSE, has also naturally evoked strong interest in soli- 
tons and their dynamics in cigar-shaped BEC [HI, |24| . 
Dark solitons in the repulsive regime and bright soliton 
and soliton trains in the attractive sector have been ex- 
perimentally produced [l5l - fl8l |. Recently, the grey soli- 
ton, having complex profile has been produced through 
soliton collision in BEC [2l| . Grey solitons have been the 
subject of some recent studies [25, 26]. These are analogs 
of complex envelope Bloch solitons in condensed matter 
systems, which can connect differently ordered domains, 
without passing through the normal phase. This natu- 
rally makes it energetically favorable, as compared to the 



Neel soliton, wherein the soliton profile, being of the hy- 
perbolic tangent type, passes through the normal phase, 
at the vanishing point of the above profile. It is interest- 
ing to note that the speed of the soliton can vary from 
zero to a maximum, depending on the depth of the grey 
soliton. Furthermore, the energy shows a maximum value 
as a function of momentum. 

The doping of nonlinear fibers with suitable multilevel 
systems has opened up the exciting possibility of natu- 
rally combining the the atomic system with NLSE, which 
can possibly lead to effective soliton control and manipu- 
lation. Recently, the localized solitons of a nonlinear fiber 
with a three level dopant has been investigated, wherein 
exact solitons of tanh — sech type have been identified 
[27I - I29I ] . The constraints of this coupled nonlinear sys- 
tem has led to strong restrictions on the pulse dynamics. 
It was found that, the frequencies have to be perfectly 
matched for these solutions to exist [l|. In this paper, we 
present exact grey soliton profiles in this system, which 
removes the above restriction on the pulse propagation. 
The complex profile's phase degree of freedom plays a 
crucial role in relaxing the constraint of frequency match- 
ing. 

The paper is organized as follows: In the following, 
grey solitons, with complex envelope, are obtained as ex- 
act solutions of the coupled nonlinear system describing 
doped optical nonlinear fibers. It is explicitly shown how 
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FIG. 1: A three level A atomic system whose transitions 1 -H> 2 
and 1 -O- 3 are in near resonance with the two modes of the 
optical fiber. The Rabi frequencies of the two modes are 2G 
and 2g. A denotes the single photon detuning. 
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FIG. 2: Stable propagation of (a) bright and (b) dark solitons 
in a three level medium in presence of nonresonant nonlinear- 
ity and fiber dispersion. This is in the regime of opposite signs 
of the group velocity dispersion at the two mode frequencies. 



the phase degree of freedom removes the restrictive equal 
frequency conditions on the tanh — sech paired pulses 
The cross-phase modulation is found to play a crucial 
role, in controlling the speed of the grey soliton. For 
the two level system, it is observed that grey soliton is 
strictly forbidden and only the SIT soliton can propagate 
in this dynamical system. Subsequently, we present the 
numerical techniques and the results of numerical prop- 
agation of our solution, which clearly shows stability of 
the obtained solutions in a wide parameter range. We 
then conclude with directions for future work. 



II. COUPLED GREY SOLITON SOLUTIONS 

We investigate the possibility of generalization of the 
solution obtained in [lj for shape-preserving pulses in a 
Kerr nonlinear two-mode optical fiber doped with 3-level 
A atoms as shown in Fig([T|). For the earlier obtained 
solutions, the frequency shift was assumed identical for 
the two possible fiber modes. The modes of the fiber are 
near resonant with the transitions of the atomic system. 
The two modes of the optical fiber may be described by 
the profile, 

Ei(z,t) = Aie-^-^ +c.c.(i= 1,2) (1) 

Ai being the slowly varying envelope, u>i and ki, the 
carrier frequency and wave number respectively of the i th 
mode. The time variation of the population of atomic 
levels are governed by optical Bloch equations. If Q(z, t) 
denotes the probability amplitude of the i th atomic level 



(Fig 1), then the following equations are obeyed within 
the rotating wave approximation 

Ci = -iA& + iGQ 2 + tgC 3 (2) 

Ca = iG^i (3) 

6 = ig'Xi (4) 

The Rabi frequencies 2g and 2G for the two field modes 
are related to the slowly varying amplitudes Ai, accord- 
ing to the relations 



where dij is the transition dipole moment matrix ele- 
ment. We take into account the effects of Kerr nonlinear- 
ity (proportional to the square of the electric field) and 
dispersion, in order to truly describe the spatio-temporal 
variation of the light pulses through the fiber. Using the 
approximation of slowly varying envelope one can cast 
the nonlinear Schrodinger equation, in Rabi frequencies, 
as 




FIG. 3: Growth of instability for (a) bright and (b) dark 
solitons in a three level medium in presence of nonresonant 
nonlinearity and fiber dispersion. This is in the regime of 
same signs of the group velocity dispersion at the two mode 
frequencies. 
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where the terms with /3j are due to group velocity disper- 
sion (GVD), those with 7, and rji represent the effect of 
Kerr nonlinearity, and the coupling of the i th mode with 
the atomic system, respectively. A dark-bright coupled 
soliton pair was reported in [lj in this system. There sta- 
ble propagation of bright and dark solitons was reported, 
in the regime of opposite signs of the group velocity dis- 
persion at the two mode frequencies. Correspondingly, 
unstable propagation was reported, when the GVDs were 
of the same sign. Figs.® and ([3]) show the solutions ob- 
tained by them in stable and unstable regimes of GVD. 



III. ANSATZ 

Motivated by the observed grey solitons in BEC 
[23l |25| , we assume the following ansatz for the complex 
envelope solitons 



g = a cos(0)sech 
G = b 
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e i(pi-p 2 )z-i(n 1 -n 2 )t 
(7d) 

(1 - a)T (7e) 



U gives the envelope velocity in the moving frame, a 
appearing on the RHS of (|7a|) gives the temporal width 
of the pulse. We have removed the restriction Oj = £^2- 
The subsequent sections show that we have been able to 
derive analogous expressions as in [l| for the unknown 
variables introduced in the ansatz. 



IV. CONSISTENCY CONDITIONS 

The ansatz is put into the Bloch and the nonlinear 
Schrodinger equations (2-6) and the coefficients of sech, 
tanh, sechtanh, etc. are collected on both sides to yield 
the following relations. The Bloch equations (2-4) give 



2 .2 
a — b 



sin(0) = a(fti - n 2 



i 2 a 2 + i(A - fii + b 2 asm(9))(a/T) 



(8a) 
(8b) 

(8c) 



The probabilities of occupation of the three levels should 
add up to one. Imposing this condition we get, 
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The nonlinear Schrodinger equations (5-6) yield 
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Equation (jlOdft can be alternately cast as, 



u = 2f3 2 n 2 + (b 2 - 2a 2 ) l2 asm(9) + r) 2 -^r (H) 

For stable propagation of pulses /?i ^ (3 2 lj. The follow- 
ing constraint can be obtained from equations (|10c|) and 



t-rim - m) + 2«i(/3i - &) = o or n x - 



(12) 

Equations (|8a|). (|T0el) and (flOfj) lead to ([131), which is an 
important result that needs to be satisfied by the non- 
linearity and dispersion in the medium at the two mode 
frequencies: 



3 + 2^ + 2^ 
7i 72 







(13) 



Another important point to note is that, for = 0, 
in order to have a solution, f2i = f2 2 is required. This 
is because of the consistency relation (|8b[) viz., sin0 = 
<r(f2i — ft 2 ). If ^1 7^ ^2) then a — and the arguments 
of all the sech and tanh functions in the ansatz (I7all7el) 



blow up to give constant solutions. Therefore Qi = Vl 2 
is required when = 0. Thus in [l|, while looking for 
pure dark soliton solution for G, the condition fii = SI2 
was automatically enforced. In comparison, in our so- 
lution, we fix Q, 2 to any value and find the rest of the 
parameters in terms of 9 (9 ^ 0, fli ^ il 2 ). There- 
fore our ansatz (|7all7e|) represents valid solutions for all 
9, given the consistency conditions arc satisfied, showing 
that grey solitons of varying depths can be solutions of 
the system. 
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FIG. 4: Grey solitons are usually velocity restricted. Often 
the relationship is as simple as it oc sin 8. In our case it is 
given by (jTTJ . The parameters used are /3i = 1, $2 = —2.5, 
71 = 1, 72 = 1, Vi = lj Vz = 1-2 and A = 1. 



The velocity of the solitons (u) has a nontrivial depen- 
dence on 9. This physically means that the velocity of 
the solitons is dictated by the depth of the grey soliton in 
the system. Fig (UJ) shows that the solitons are velocity 



restricted with maximum value attained at 9 = tt/2. 



V. STABILITY OF THE COUPLED GREY AND 
BRIGHT SOLITON PAIR 

A. Linear stability analysis 

Linear stability analysis is done by perturbing around 
the known solutions. Perturbing around 5 complex wave- 
forms results in 10 coupled differential equations in 10 
variables with complicated potentials (composed from 
the solutions themselves), which are difficult to handle 
either analytically or numerically. Hence, we have car- 
ried out numerical evolution of the pulse for the pur- 
pose of stability check. Here, we state the equations of 
motion of the perturbations, as derived from the equa- 
tions of motion for g, G and Q (2-6) by substituting 
#->• (g + 5g)e l ( piz - nit \ G ->■ (G + SG)e l(p2Z ~ n2t \ Ci -> 
(Ci+£C~i)e i(pi *" nit) , C2 -> (C2+SC2)e i(pi - p2 >e- i( - Ul - u ^ t 
and (3 — > ((^3 + 5^3) . The quantities g, G and Q are given 
by the RHS of equations (|7all7ep respectively, without the 
exponential terms. The perturbations satisfy, 
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where the operators Tij stand for: 



B. Numerical stability analysis 



F33 
T44 

Tl2 

r 2 i 



- -i(d z + ipi) + /3i(<9 t 2 - 2itt 1 d t - 9. 

27i(M 2 + |G| 2 )- 7 iS 2 conj 
■ -i(d z + ip 2 ) + /3 2 (5 2 -2ifl 2 d t 

27 2 (M 2 + |G| 2 )- 7 iG 2 conj 
: d t - ifl + iA 

■■d t -i(ni-n 2 ) 

■ -2-fxgG* - 2715G conj 
: -2-f 2 g*G - 2j 2 gG conj 



(14a) 

(14b) 
(14c) 
(14d) 
(14e) 
(14f) 



Here, conj is the function that takes the complex con- 
jugate of the term it acts on eg., conj(<5(jj) = Sg*. The 
above matrix equation is thus a condensed form of an 
actual set of 10 coupled partial differential equations in 
10 variables. 



In this section, we investigate the propagation dynam- 
ics of the grey and bright solitons by integrating the full 
set of coupled nonlinear Schodinger and Bloch equations. 
We use Runge-Kutta and split step operator methods to 
investigate the spatio-temporal evolution of the optical 
pulses in fiber doped with three level A system. The 
initial shapes of the input pulse and the initial atomic 
populations can be obtained in the limit of z — and 
t = from the analytical expressions (|7all7b[) and (|7cll7c[) 
respectively. We use grey and bright solitons with width 
3.6712 psec in the present investigation. Parameters used 
in the numerical simulation obey self consistent relations 
(fT2"| and (fP3")l . We first explore the stable propagation 
of the grey and bright solitons in presence of GVD with 
same character. Figure <j5j) shows that the amplitudes of 
the output pulses are retaining their initial shape. Both 
solitons propagate through the medium with same group 
velocity. A small oscillation though has been noticed at 
both the ends of the grey soliton due to the modula- 
tion instability. Our analytic and numerical solution of 
grey-bright soliton pairs are presented in the inset of the 
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FIG. 6: Amplitude and phase of the dark and grey solitons 
are plotted as a function of time at the exit of the medium. 
All other parametrs are as in Fig. [S] 



dark soliton develops instability, both in intensity as well 
as in the phase while grey soliton propagates through the 
medium without loss of generality upto sufficient length 
of propagation. Hence the instability due to the same 
GVD sign can be compensated with the proper choice of 
the blackness parameter 6. 

Figures (|3i)-([7p) show the results for the spatio- 



FIG. 5: Spatio-temporal evolution of the grey-bright soli- 
ton pair in medium plotted against the spatial length of the 
medium at different propagation times. The other parameters 
used are A=0.33361, B=0.19261, pi = 1.0039, p 2 = 0.037098, 
Qi = 0.15584, 0,2 = 0,9 = 0.60908, /3i = -0.25, f3 2 = -1.25, 
71 = 1, 72 = 1, A = 1. u = 1.4805, a = 3.6712, V = 1.0338, 
rji = 1, and 772 = 1.2. 



figures |[SJ) at the propagation length z = 3m. The nu- 
merical results obtained from the Runge-Kutta and split 
step operator methods match well with those obtained 
analytically. Next we study how the blackness parame- 
ter 9 leads to stable propagation of soliton through the 
medium even though GVD has same sign. A comparative 
study of the dark and grey solitons through the medium 
is shown in the figure |U The dark soliton can be ob- 
tained with 9 = whereas < 9 < tt/2 represents the 
grey soliton. For dark soliton, the intensity of the dip 
always vanishes but the same for the grey soliton always 
remains. The variation of phase with time with propa- 
gation length for both dark and grey solitons is seen to 
be antisymmetric in nature. Interestingly, dark soliton 
shows a sudden change in phase across the line center 
while grey soliton displays no change in shape of phase 
across the line center during the length of propagation. 
In the course of propagation, the antisymmetric phase 
has a temporal shift due to the presence of self and cross 
phase modulations of the medium in both the cases. The 




FIG. 7: Spatio-temporal evolution of the grey-bright soli- 
tons using parameters such that the group velocity dispersion 
have the opposite sign. We used A=0. 29551, B=0.25592,pi = 
0.77581, p 2 = 0.065494,^1 = 0.072505, Q 2 = 0, 9 = 0.51291, 
/3i = 1.0, (3 2 = -2.5 ,71 = 1, 72 = 1, A = 1. u = 2.6827, 
a = 6.768, T = 1.0337, 771 = 1, and 772 = 1.2. 
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temporal evolution of grey and bright solitons in presence 
of the opposite group velocity dispersion. Both the grey 
and bright solitons maintain their initial shape. 

Conclusions 

In conclusion, grey solitons, analog of Bloch solitons in 
condensed matter systems, have been identified as exact 
solutions of the coupled nonlinear Schrodinger and Bloch 
equations. The complex nature of the profile played a 
crucial role in removing the restrictive frequency match- 



ing condition of tanh— sech paired solutions We nu- 
merically evolved our solution and found it to be stable in 
both regimes of group velocity dispersion. We note that, 
for the case of coupling with a two-level system (g — > 0, 
Ca — * 0) instead of the A system considered here, the 
pure tanhp and grey solitonic profiles are prohibited for 
G, while a pure sech profile is allowed. The fact that our 
solutions are stable for some parameters in both regimes 
of group velocity dispersion signs, is very interesting and 
needs further investigation. 
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